


♦ 









I 





NAVAL POSTGRADUATE SCHOOL 

Monterey, California 








A NUMERICAL STUDY OF RAIN-INDUCED 
SURFACE GRAVITY WAVE ATTENUATION 




by 




David W. Howell 
• » <» 




June 1989 


Thesis Advisor: 


Jeffrey A. Nystuen 



Approved for public release; distribution is unlimited. 



T9iiADQ7 









■I 




t 





.'ndassified 

scunty classification of this page 



REPORT DOCUMENTATION PAGE 



la Report Security Classification Unclassified 



Ib Resirictise Markings 



2a Securiiv Classification Authority 



2b Declassification Downsradins Schedule 



3 Distribution Availability of Report 

Approved for public release; dislribulion is unlimited. 



J Performing Organization Report Numberfsi 



5 Monitoring Organization Report Number(s) 



Sa Name of Performing Organization 

Naval Post 2 Taduate School 



6b Office Symbol 
(if applicable) 35 



7a Name of Monitoring Organization 

Naval Postgraduate School 



be Address (cl:y, state, and ZIP code) 

Monterev. CA 93943-5000 



7b Address (city, state, and ZIP code) 

Monterev' CA 93943-5000 



sa Name of Funding Sponsoring Organization 



Sb Office Symbol 
{ if applicable ) 



9 Procurement Instrument Identification Number 



>c Address (dry. state, and ZIP code) 



10 Source of Funding Numbers 



Program Element No Project No Task No Work Unit Accession No 



11 Title { include security classification) A NUMERICAL STUDY OF RAIN-INDUCED SURFACE GRAVITY WAVE AT - 
FENUATION 



2 Personal Author(s) David W . I 


owell 


3 a T ype of Report 


:3b Time Co\ercd 


14 Date of Report {year, month, day) 


1 5 Pace Count 


Master's Thesis 


From To 


June 19S9 


3S 



6 Supplementary Notation The views expressed in this thesis are those of the author and do not reflect the official policy or po- 
ition of the Department of Defense or the I'.S. Government. 



7 Cosati Codes 



"lejc 


G.^-oup 


Subgroup 















IS Subject 1 erms ( continue on reverse ij necessary and identify by block number) 

Ocean Waves, Attenuation, Vortex Rines 



9 Abstract {continue on reverse ij necessary and ideniijy by block number) 

Strong rain-induced mixing in a thin surface layer is numerically shown to greatly increase surface gravity wave atten- 
jation. This case study uses a single wavelength (2.S m) together with two mixed layer depths (10 and 20 cm). The rain- 
nduced mixing is simulated by var>ing kinematic viscosity within the mixed layer from 10”Mo . molecular to 

•trong turbulent mixing, respectively. The results indicate that surface gravit\* wave attenuation in the presence of a thin 
ain-induced mixed layer can increase by a factor of up to 6000 times the attenuation rate due to molecular viscosity alone. 
I’his indicates that rain need onl\ mix the top 10-20 cm surface la\ er to effectively dampen short surface gravity waves. 



0 Distribution .Availability of Abstraci 

N unclassified unlimited □ same as report D DTIC users 



21 Abstraci Security Classificaiion 

Unclassified 



la Name of Responsible Individual 

■elT Nvstuen 



22b Telephone {include Area code) 

(40S) M6-2917 2552 



22c Office Symbol 



D 1 OR\I 1473.S-: MAR 



S3 .APR edition may be used until exhausted 
All other editions are obsolete 



security classification of this page 

Unclassified 



1 



Approved for public release; distribution is unlimited. 

A Numerical Study of Rain-Induced 
Surface Gravity Wave Attenuation 

by 

David W. Howell 
Lieutenant, United States Navy 
B.S., Texas A&M University, 1979 

Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN METEOROLOGY AND PHYSICAL 

OCEANOGRAPHY 

from the 



NAVAL POSTGRADUATE SCHOOL 
.lunc 1989 



TABLE OF CONTENTS 



I. ABSTR.ACT 1 

II. INTRODUCTION 2 

A. SIGNIFICANCE 3 

B. LITEIUATI RE REVIEW 5 

E Surface Effects Due to Rain 5 

2. Characteristics of Vortex Rings 6 

3. Rain-induced Mixed Layer Changes 7 

4. Manton's Analytical Approach S 

III. .METHODS 10 

A. THE .MODEI 10 

B. THE WAVE 14 

C. THE .MIXED LAVER FROBLE.M 22 

1. Single Layer Wave 24 

2. Rain's Elfects on Wa^cs 25 

a. lu cm Mixing Depth 25 

b. 2'i cm .Mixing Depth 25 

IV. CONCLUSION 26 

REFERENCES 2S 

INITIAL DISTRIBUTION LIST 30 



LIST OF TABLES 



Table 1. RESULTS FOR A 3 CM AMPLITUDE WAVE 21 

Table 2. RESULTS FOR A 6 CM AMPLITUDE WAVE 21 

Table 3. RESULTS FOR A 3 CM AMPLITUDE WAVE, 10 CM MIXING 

DEPTH 24 

Table 4. RESULTS FOR A 6 CM AMPLITUDE WAVE, 10 CM MIXING 

DEPTH 24 

Table 5. RESULTS FOR A 6 CM AMPLITUDE WAVE, 20 CM MIXING 

DEPTH 24 



iv 



LIST OF FIGURES 



Figure 1. Mixed layer depth prediction 4 

Figure 2. Typical cell configuration 11 

Figure 3. Kinematic viscosity grid configuration 13 

Figure 4. Initial wave configuration for all experiments 15 

Figure 5. Wave propagation through the tank 16 

Figure 6. Wave propagation -■ 17 

Figure 7. Natural log of u vs depth 18 

Figure 8. Potential and Kinetic energy vs time 20 

Figure 9. Potential and Kinetic energy vs time 22 

Figure 10. Rate of energy change 2 



V 






ACKNOWLEDGEMENTS 



1 would like to thank Professor Jeffrey Xystuen and Professor Roland Garwood for 
their unending support and answers to my many questions. A special thanks also goes 
to my wife, Anna, for her support, guidance, skillful typing and patience in the writing 
of this manuscript. 



VI 



I. ABSTRACT 



Strong rain-induced mixing in a thin surface layer is numerically shown to greatly 
increase surface gravity wave attenuation. This case study uses a single wavelength (2.8 
m) together with two mixed layer depths (10 and 20 cm). The rain-induced mixing is 
simulated by varying kinematic viscosity within the mixed layer from 10”‘ to 10’^ m- s~' , 
molecular to strong turbulent mixing, respectively. The results indicate that surface 
gravity wave attenuation in the presence of a thin rain-induced mixed layer can increase 
by a factor of up to 6000 times the attenuation rate due to molecular viscosity alone. 
This indicates that rain need only mix the top 10-20 cm surface layer to effectively 
dampen short surface gravity waves. 



II. INTRODUCTION 



For centuries sailors and casual observers living along the coast have observed that 
ocean waves are dampened, and the sea becomes flat during rain. This dampening is 
evident especially during rain events when winds are light or calm. This dampening must 
be a result of energy transfer from a raindrop that strikes the ocean surface. The 
raindrop impact may generate capillarv' waves, subsurface turbulent mixing, and sound 
or heat. Only subsurface turbulent mixing is likely to interact with surface gravity 
waves. This process, first suggested by Reynolds (1900), takes the form of vortex rings 
which propagate downward, mixing the water in a near surface layer as the rings disin- 
tegrate; thus enabling the vortex rings to e-Hcctively increase the turbulent mixing in a 
thin layer near the surface of the water. 

Manton (19' 3) used an analytical approach based on Reynolds' idea, to argue that 
rain-induced mixing will cause surface gravity wave attenuation. Manton assumed a 
constant eddy viscosity throughout the water column to solve his governing momentum 
equations, 'fhe assumption of no depth dependence of the mixing seems unrealistic. 
The scope of this study is to examine, by numerical methods, the efi'ects of increased 
turbulent mixing due to rain in a thin surface layer. It is hypothesized that by increasing 
the turbulent mixing in a thin surface layer, significant attenuation of a surface gravity 
wave will result. This study focuses on waves having a wavelength of 2.8 m and ampli- 
tudes of 3 and 6 cm. Eddy viscosity is varied from 10"‘ to lO’Mn ^ . The mixed layer 

depth is either 10 cm or 20 cm. The rate of kinetic energy change at each time step was 
calculated and used to determine the rate of wave attenuation. The use of a numerical 
calculation allowed viscosity to be varied with depth. Alternatively. Manton employed 



a velocity potential solution; however, this precluded solving for the rotational Hows 
occurring here. 

A. SIGNIFICANCE 

The dampening of waves by rain can influence satellite remote sensing observations 
at the sea surface, in addition to changing the characteristics of the water column. Sat- 
ellite borne scatlerometers and Synthetic-Aperture Radars (SAR) measure microwave 
energy backscattered to space from the ocean's surface. Scatlerometers are able to de- 
termine wind speeds based on the energy backscattered from the sea surface by short 
surface gravity waves (appro.ximately 30 cm). If these waves arc attenuated by rain. 
scatterometer-deri\ ed winds will be in error for areas under a rain event, giving the im- 
pression that there are light winds present when in fact greater winds may be present. 

Fu and Holt (19S2) noticed that squall storm lines were easily detected in S.-\R im- 
ages. They found that down draft and rain areas were delineated by changes in image 
brightness. The image was brighter in the outflow area due to a rough sea surface. A 
smooth dark image caused by the absence of backscattered energy suggested that in 
areas of rain the capillary and short gravity waves had been dampened. 

Rain also innuences the water column. Livezey (19S8) found that rain events could 
cause changes to the seasonal mixed layer depth indirectly effecting heat loss to the at- 
mosphere. Using a total kinetic energy budget model. Garwood (1989) demonstrated 
that a rain-generated iruxed layer would form in a stratified water column. .A 2t)0 cm 
mixed layer depth (.MUD) was produced when turbulent kinetic energy (TKE) of a rain 
event was injected into the upper layer of the ocean, see Fig. 1. Within 1000 seconds 
the .MLD is seen to be approaching steady state. A continued deepening of the .MLD 
will occur, but at a much slower rate until equilibrium has been reached. The results are 
tentative, but do suggest that rain has a possible role in ocean mixed layer dynamics. 
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1 line (seconds) 

Figure 1. Mixed layer depth prediction: A 200 cm mixed layer has formed by 

rain-injected TKE. 

When a raindrop strikes a fluid surface there are other interesting elTects. One effect 
is a roughening of the sea surface from depression cusps that break the smoothness of 
the sea surface. In conjunction with the wind, the rougliened surface caused an increase 
in the surface stress. This stress increase augments the transfer of energy into the w ater 
column. Caldw'cll and Elliott (1970) found that a rain rate of several centimeters per 
hour in light winds could produce stresses comparable in magnitude to the wind stress 
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for winds as high as 10-20 ms-‘ without rain. A hea\y rain can also produce bubble 
clouds, a layer of bubbles beneath the wave surface, similar to those produced by 
breaking waves. Thorpe (1984) used these bubble clouds to examine the effects of 
Langmuir circulation by following their progress through the water column. In his study 
he found that the Langmuir circulation increased the flux of gases down into the water 
column. Rain is thus able to effect directly the surface roughness and stress, and to ef- 
fect indirectly the flux of gas into water column. 

B. LITERATURE REVIEW 

1. Surface Effects Due to Rain 

Reynolds (1900) observed water drops falling into a glass vessel and he postu- 
lated the first explanation for wave attenuation. From his study of vortex rings, he be- 
lieved that the *’ knocking down ’’ of waves resulted from mixing of subsurface water by 
vortex rings formed by the raindrops hitting the surface. As the vortex ring descended, 
it entrained fluid from the moving surface layer and carried this fluid downward. 
Reynolds reported that the nonmoving fluid from below would rise up and replace the 
surface fluid, diminishing the motion of the surface which is essential to the continuity 
of the wave. 

Later. Ashton and O'Sullivan (1949). Barnaby (1949), and Sainsbun.' and 
Cheeseman (1950) further investigated Reynolds' theory. Each of these groups con- 
ducted experiments involving drops of dye striking a water surface. Their conclusions 
v’re inconclusive in testing Reynolds' theory. Although vortex rings were observed in 
all experiments, the rings' effect on wave attenuation was not clearly established. 

The interaction of water drops on surface waves was explored subsequently by 
Siscoe and Levin (1971). They found two ven.' different processes, depending upon the 
raindrop's point of incidence on the wave; a reflection process and an absorption proc- 
ess. In an absorption event, the drop strikes at or near the wave crest and TKE is in- 



corporatcd into the wave field. A reflection event occurs in the event of rain striking the 
trough, and the TKE of the drop goes into producing subsurface turbulent mixing. They 
found that during an absorption event the energy supplied by an impacting drop ampli- 
fied the surface wave field which propagated out from the point of impact. Drops were 
colored by dye to investigate further wave formation and subsurface interaction. The 
dye did not penetrate the surface for an absorption event, while during reflection events, 
dye penetrated to the subsurface waters due to drop-induced mixing. The observed 
motions of the dye that accounted for this mass transfer were vortex rings and weakly 
defined circulation cells that reached a depth of 5-10 cm with an occasional ring pene- 
trating deeper. 

2. Characteristics of \'ortex Rings 

When raindrops strike the sea surface a portion of their energy is transferred 
into the water column in the form of vortex rings. Kecdy (1967) conducted an extensive 
study on the formation of vortex rings by free surface interactions. He allowed water 
drops of various radii to fall from heights of 2.5 inches or less, and he observed ring 
formation and penetration. Kcedy found that vortex rings penetrate to a maximum 
depth of approximately 10 cm. Also noted were density differences between the drop 
and the base solution which were important in determining the maximum penetration 
depth of the vortex ring. For a drop with a density lower than the density of the base 
solution, the vortex ring expanded in diameter more rapidly than if both drop and base 
solutions were of the same density. 

-Vlaxworthy (1972) observed the propagation of vortex rings produced using a 
membrane apparatus in an experimental tank. He found that as the rings traveled they 
continuously entrained surrounding fluid and grew in volume. At the same time, the 
rings were losing vorticity into a wake and decreasing in velocity. Maxworthy also ob- 
served two different processes when two similar rings, one ahead of the other, came into 
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contact. If both rings had the same velocity, the rearward one would distort and squeeze 
into the forward ring through the rear wake area. The rearward ring was then wrapped 
around the inner part of the forward ring to form a single larger vortex. If the velocity 
of the rearward ring was larger than that of the forward ring, the rearw'ard ring's entry 
into the forward ring caused the combined ring to oscillate longitudinally. If the trailing 
ring's velocity is much greater than that of the front ring Maxworthy found that the 
newly formed ring would disintegrate with a faster ring propagating ahead of a slower 
moving ring. 

Batchelor (1967) had previously found a similar occurrence involving two similar 
vortex rings and suggested an interaction analogous to “ leap frog. ” The velocity field 
of the rear ring has a radially outward component at the position of the forward ring. 
As the radius of the front ring increases, travel speed decreases. At the same time, the 
rearward ring undergoes a relative increase in speed and passes through the larger 
vortex, becoming the lead ring. This maneuver is then repeated. The interaction of 
these rings suggests that random rings produced in nature will collide, changing the 
ring s characteristics and velocity. The gross effect will be to reduce the penetration 
depth expected from a single vortex ring. 

3. Rain-induced Mixed Layer Changes 

Rather than study individual drops and vortex rings. Green and Houk (1979) 
studied populations of raindrops and the subsequent mixed layer formation. Observa- 
tions of changes to the mixed layer depth (.MLD) were made by determining the location 
of the thermoclinc. This circumvented the need to consider individual vortices in favor 
of an integrated turbulent mixed layer. Green and Houk found that a mixed layer had 
developed after 20 minutes of rain, with a deeper MLD associated with the larger drop 
sizes. The depth of the mixed layer in salt water was approximately one third that of 
freshwater due to the larger downward buoyancy flux. All other parameters were held 



constant. Green and Houk determined that the majority of the mechanical energy flux 
went directly into subsurface turbulence. The overall effect was a deepening of the MLD 
as the energy from raindrops diffused down into the water column, deepening the MLD 
10-30 cm in freshwater and 5 to 10 cm in saltwater. 

4. Manton's Analjlical Approach 

Manton (1973) presented an analytical approach to explain rain-induced vortex ring 
mixing and the resulting wave attenuation. Through a scale analysis, using the 
raindrop s terminal velocity (V) and the horizontal wind speed (U), Manton defined a 
vertical and horizontal eddy viscosity: 



v,=/M7 (2.1) 

( 2 . 2 ) 



where/ represents the volume fraction of drops in the air; and / represents the turbulent 
length scale determined by the average distance between drops in the air. In order to 
use vertical and horizontal eddy viscosities in the momentum equations, Manton needed 
to prove that equations 2.1 and 2.2 were independent of depth. Through analytical 
means he proved that the shear layer thickness for a wave, on the order "'as 

smaller than the penetration depth of vortices (10-20 cm). Erroneously, the shear layer 
thickness is the depth associated with friction from a boundary and not the depth asso- 
ciated with the wave energy, 1 k. By assuming that the maximum mixing depth is related 
to vortex ring penetration, mixing should only occur at a depth of 10 to 20 cm, which 
is much shallower than the depth associated with the wave energy. Manton solved the 
linearized momentum equations using a potential velocity solution, implicitly assuming 
irrotational flow. However, when mixing is present the flow becomes rotational. It is 
still possible to predict the viscous decay of the waves by assuming that the mixing is 



uniform over depih (Lighihill, 1978). If and v,, are assumed to be equal, .Manton's 
results would agree with such theory; however, Manton used separate values of 
v'v and V;, to incorporate the influence of wind into his solution. This is an inappropriate 
use of the viscous decay terms in the momentum equations. For a vorte.x ring, the 
horizontal and vertical mixing scales are of the same order, thus and v* should be re- 
placed by a single eddy viscosity coefficient, v , for the mixed layer. Wind, which is not 
considered in this study, is normally introduced as an additional surface stress source 
term in the momentum equations and does not appear as part of the viscous dissipation 

term. Thus there are two serious objections to .Manton's effort; 

1. introducing mixing over the entire water column, and 

2. incorporating the influence of wind into v. 

This study uses a numerical solution of the Navier-Stokes equations for a surface gravity 
wave. A single viscous coeflicient v will be applied to a thin " mixed layer " and the 
resulting attenuation will be estimated and compared to classical theory (Lighthill. 197S). 
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III. METHODS 



A. THE MODEL 

A numerical model, SOLA-VOF (Nichols et al, 1980), developed at the Los Alamos 
National Laboratory' in New Mexico was used for this study. This model is a general 
two-dimensional finite differencing Eulerian code which allows for complicated compu- 
tational domains with variable boundary conditions, free surface shapes, variable surface 
tension and compressibility. Behavior of the fluid is determined by solving the Navier- 
Stokes equations: 
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where u is the horizontal velocity, v is the vertical velocity, g is the body acceleration due 
to gravity, v is the coefficient of kinematic viscosity, and p is the fluid density. In addi- 
tion to the momentum equations the incompressible mass conservation equation: 



4 ^ + 4 ^ = 0 

cx cy 



(3.3) 



must be satisfied at each time step. 1 o advance a solution through a one time increment 

the model must go through the following three steps: 

1. Explicit approximations of Eq. (3-1) and Eq. (3-2) are used to compute the first 
guess for the new time-level velocities using the initial conditions or the previous 
time-level values for all advective, pressure, and viscous accelerations. 

2. To sati'^fy the continuity equation. Eq. (3-3), pressures arc iteratively adjusted in 
each cell and the velocity changes induced by each pressure change are added to 
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the velocities computed in step (1). A repetitive process is required because tlie 
change in pressure needed in one cell to satisfy Eq. (3-3) will upset the balance in 
the four adjacent cells. 

3. The cell fluid function defining fluid regions is then updated to give the new fluid 
configuration. 

In this manner the lluid solution can be advanced through all desired time intervals. 
Horizontal and vertical velocity (m,v), pressure (P), and fractional volume of fluid in a 
cell (F) are calculated at each time step at the cell positions shown in Fig. 2. 




Figure 2. Typical cell configuration: U and V velocities are found on the sides 

and top of the cell with pressure and fluid percentage of each cell found 
in the center of the cell. 
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The fractional volume of fluid function (F) is used to identify grid cells that contain 
fluid, and to delineate surface, interior or boundary cells. A free surface cell is defined 
as a cell containing a nonzero value for F and having at least one neighboring cell that 
contain a zero value for F. Cells with nonzero F values and no empty neighbors are 
treated as cells full of fluid. 

Introduction of variable viscosity into the code required a change to the code. This 
was accomplished by; (1) establishing a viscosity grid and (2) encoding this viscosity grid 
into the viscous terms. The viscosity grid was established by identifying three areas, v,, 
the mixing layer at the surface of the wave, Vj, a transition zone from the mixed layer to 

Vj ”j” Vj 

the turbulent free layer defined as — ^ ^' 3 ’ layer, (effected by molecular 

viscosity only). The turbulent mixing layer was defined as those cells which were be- 
tween the free surface and the desired mixing depth. An example is given in Fig. .3 de- 
picting the mixing layer following the surface. The height change is from the numerical 
initialization of the wave. 



J 

V’l 



Figure 3. Kinematic viscosity grid configuration: The setup of the viscosity grid 

where v, is the surface mixing layer, vy, the transition layer, and v, the 
base fluid layer. 



The viscosity grid was redefined at each time step resulting in a moving mixing layer 
having a constant thickness. A four corner average was then used to obtain a value of 
V for each computational cell. The single derivative form of the momentum equations 
(Harlow 1980): 
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was used to reduce numerical noise associated with calculating second derivatives at the 
boundary’ of the two layers. 

B. THE WAVE 

This study focused on the attenuation of deep water waves resulting from rain- 
induced mixing. To simulate this a computational domain was initialized representing 
a fluid vessel or tank 1.4 m wide and 3.0 m deep. A Mean Fluid Height (MFH) of 2.8 
m was selected. The tank depth of one wavelength was chosen to ensure that bottom 
effects would be negligible during the propagation of the wave. According to deep water 
wave theory, velocity at the base of the tank is smaller than the surface velocities by a 
factor of c"*' = 1 . Furthermore, a small amplitude wave {3 and 6 cm) was chosen 

to maintain stability of the wave. During the early part of the project the wave was 
found to break if the amplitude became greater than 10 cm due to the manner in which 
the code assigned velocities at the free surface. Normally a surface velocity is set to be 
the same value as the nearest horizontal neighbor cell. Depending on the direction of 
wave propagation this horizontal movement of velocities occasionally produced an un- 
realistic velocity discontinuity in the vertical. This problem was corrected by making the 
nearest neighbor cell directly below the surface cell. 

.‘\ 2.S m wavelength standing wave was chosen to minimize continuation boundary 
problems associated with a progressive wave, and to facilitate a simple initial condition 
having no velocities. A standing wave is actually two progressive waves propagating in 
opposite directions given by: >/ = cos kx cos cor where 7 /q is the wave amplitude, k is 
the wave number, and oj is the radial frequency. A progressive wave was simulated by 
initializing one-half wavelength of the standing wave using the instantaneous wave 
height: ?; = MFH — ij^coskx at t = 0 for a domain of 0 < x < -^. The waves initial 
jnfiguration is seen in Fig. 4 under steady state conditions. Velocity vectors are also 
plotted for the initial steady state condition and are represented by a vector of zero 
length. 
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Fi gure 4. Inilial Ma\ e conligui ation for all experiments: I hc one-half wavelength 

at rest prior to being released. 



The wave was released from rest and observed as it propagated through the tank. Wave 
dynamics were reviewed and compared to theoretical values to ensure that the wave and 
code had been properly developed. 

Validation of the wave was accomplished by examining the period, the vertical U 
velocity profile, the potential and kinetic energy, and the attenuation characteristics of 
the wave. According to deep water wave theory a period of 1.34 s is derived. In 
Tig. 5 and Fig. 6 the wave is seen to lake appro.ximatcly 1.3 s to propagate from a 
maximum amplitude at the right wall to the left and back to the right wall coinciding 
with theory. 
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T= 0.90 s CYCLE = 90 T= 1.00 s CYCLE = 100 T= 1.10 s CYCLE =110 





.01 rii/s 



Figure 6. 



Wave propagation: The wave is now propagating toward the rigl)t wall 

of tile tank. 
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Figure 7. Natural log of u vs deptli. 



To validate the velocity field a vertical U velocity profile was examined and compared 
to theorv': u = iig The theoretical velocity decay profile was derived and plotted, 
Fig. 7. The numerical profile was consistent with theory except for an inconsistency 
seen near the bottom of the tank which may be due to numerical noise contanrination. 

The energy of a standing wave is known to be periodic alternating from Potential 
Energy (PE) to Kinetic Energy (KE) over time. Theoretical values of the total energy 
density are obtained by adding the energy density of the two progressive waves, 
(Vj p g ijI ).) each with half the amplitude, = '/: >/o) of the standing wave. For a wave 
amplitude of, ?/o , 6 cm the total energy is V*pgu\ = 12.348 J and for a 3 cm amplitude, 
3.US7 J. To compute energy numerically, PE density was calculated from: 
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(3.6) 



PE = Vi p g rj^ dx dz 

by summing over the surface of the wave, where rj is found to be the distance the wave 
surface is above or below the .MFH. KE density was found by summing: 

KE=Vi pV'^dxdydz (3.7) 

over the domain of the tank where = (i/^ + v^) for each cell. The energy for a 6 cm 
amplitude wave propagating under normal conditions of molecular viscosity 
(v = is seen in Fig. 8. The peak energy values are consistent with theory. 

Initially PE is a maximum and KE is zero. As the wave propagates PE approaches zero 
and KE reaches a maximum at each one-quarter cycle. PE becomes a maximum at each 
one-half cycle as KE approaches zero. 
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Figure 8. Potential and Kinetic energy vs time: The sinusoidal nature of PE and 

KE is depicted for the wave as it propagates over time, v = s'h 

Kinematic viscosity was used to represent the turbulent mixing occurring in the 
water column. The values used ranged from representing molecular effects to 

10"^/?t^s"' representing high turbulent conditions. This range included values found by 
Manton: 3 x s"' for a light rain event and 1.6 x 10'^/?i^5"' for a heavw rain 

event. This range also included Green and Houk's estimate for the mixed layer eddy 
viscosity coefficient, 5 x 10"’ m’ j-'. The kinematic viscosity was kept constant over the 
entire water column to determine if attenuation effects were following theory. Both 
amplitudes were used for this test. Two difl'erent methods were used to determine the 
rate of energy change over time. The slope of the decay of peak kinetic energy was used 
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as a direct representation of the rate of energy change. From Fig. 9 the attenuation due 
to V = 10'^ m' gives a normalized slope of0.15 which is nearly consistent with theorv', 
4vk^ = 0.2. A more accurate method was found using the internal dissipation function 
(Schlichting. 1960): 




This represents the decrease of energy from viscous effects and does not include the re- 
distribution of energy into the wave. At each time step the change of energy was cal- 
culated and normalized by the instantaneous kinetic energy. Final processing of the data 
was then conducted by taking the average of the rate of change of energy divided by 
KE, i.e.,<l)/A'£. for KE > 1. This was required due to an erroneous spike developing as 
KE became less than one. as seen in Fig. 10. This was done for both the 3 and 6 cm 
amplitude cases, the results of which are given in Tables 1 and 2. showing good agree- 
ment with theor}'. 



Table 1. RESULTS FOR A 3 CM AMPLITUDE WAVE. 
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Table 2. RESULTS FOR A 6 CM AMPLITUDE WAVE. 
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.0627 


.(UMO 


l.E-2 


.384 


.4(M» 




Figure 9. Potential and Kinetic energy vs time: The sinusoidal nature of PE and 

KE are depicted for the wave as it propagates over lime, 

V = 10-2 5-'. 



C. THE MIXED LAYER PROBLEM 

Due to observations by Keedy, Siscoe and Levin, and llouk and Green two mixing 
depths, 10 anu 20 cm were used in this study. I he 10 cm depth was r presentative of 
the maximum penetration depth of vortex rings formed by raindrops striking the ocean 
surface while 20 cm represented an increased mixing depth penetration, but occurring 
more infrequently. Kinematic viscosity of the base layer was set at lO-^w^s-'. 'fhe 
surface mixing layer was varied from 10“^ to lO-^m’s"'. As for the single layer case, 
<1>/A'E was calculated at each time step and averaged for all values of KE greater than 
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Figure 10. Rale of energy change: Attenuation for <I> ( — and <\'>IKF. (---) with 

one. These results were compared to tlie expected theoretical values of the two layers. 
The theoretical values shown are where v, is the mixing if the viscosity coeflicient in 
the mixed layer is used and Vj is the mixing if the viscosity coefficient (molecular) of the 
lower layer is used. 

The effects of increasing v in the top 10 cm layer is expected to increase attenuation, 
although less than expected if the entire water column is more turbulent. Doubling the 
mixing layer to 20 cm is hypothesized to double the attenuation seen in the 10 cm sur- 
face layer case. By increasing the penetration depth more kinetic energy should be taken 
out of the wave, i.e., dampen the wave. 
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The results shown in Tables 3 and 4 are for a 10 cm mixed layer and in Table 5 for 
a 20 cm mixed layer. 



Table 3. RESULTS FOR A 3 CM AMPLITUDE WAVE, 10 CM MIXING 
DEPTH. 



Nu 


0/A'£ 


Theory’ 8v/c- v,/vj 


l.E-4 


.00162 


.0040 '.00004 


3.E-4 


.00477 


.01 20, .00004 


l.E-3 


.0157 


.0400 '.00004 


1.6E-3 


.0243 


.0640, .00004 


l.E-2 


.1269 


.4000, .00004 



Table 4. RESULTS FOR A 6 CM AMPLITUDE WAVE, 10 CM MIXING 
DEPTH. 



Nu 


(b/A'£ 


Theory Sv/<^ v,/v 3 


l.E-4 


.00165 


.0040 .00004 


3.E-4 


.00489 


.0120 .00004 


l.E-3 


.0159 


.0400, .00004 


1.6E-3 


.0252 


.0640 .00004 


l.E-2 


.1379 


.4000 .00O<:»4 



Table 5. RESULTS FOR 
DEPTH. 


A 6 CM AMPLITUDE 


WAVE, 20 CM MIXING 


Nu 


0)/A£ 


Theory SvA- v,/vj 


l.E-4 


.00259 


.0040 .o0(»04 


3.E-4 


.00773 


.0120 .00004 


l.E-3 


.0256 


.0400 .OOOi.M 


1.6E-3 


.0409 


.0640 .00004 


l.E-2 


.2440 


.4000 .00004 



1. Single Layer Wave 

Using molecular mixing, v = 10-‘ m" s"' , (Tables 1 and 2) a 3 cm amplitude wave 
required 46,000 wave periods and a 6 cm amplitude wave required 82,000 wave periods 
to reduce KE to an arbitrary cutoff value of 0.5 N-m. By comparison, the rate of energy 
change for a 3 cm amplitude wave with a 10 cm mixed layer (r, = 10^- //f 5-') showed a 
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3,000 fold increase in the rate of energy loss, requiring only 14 wave periods to dampen 
the wave. An almost identical rate of energy loss is seen for the 6 cm amplitude wave 
with dampening occurring in 22 wave periods. If the mixing depth was doubled to 20 
cm there was a 6,200 fold increase in the rate of energy change which is nearly double 
that of the 6 cm amplitude wave with a 10 cm mixing depth. Thus by increasing tur- 
bulent mixing in a thin surface layer an increase in the rate of energy change (3000-6000 
fold or greater) can be produced. This is especially noteworthy for short gravity and 
capillary waves whose energy is mostly contained within the mixing layer. 

2. Rain's Effects on \\'aves 

a. 10 cm Mixing Depth 

To model the effects of rain, three viscosities were directly associated with 

rain events: 

1. 3.0 X 10'^ m? 5"' for a light rain event (5 mmihr rain rate) (.Manton) 

2. 1.6 X 10*^ nt- s~' for a hea \7 rain event (12 mmjhr rain rate) (Manton) 

3. 5.0 X 10’'/7i-5'‘ for a hca\A rain event {lAmmlltr rain rate) (Green and Ilouk) 

Viscosities (1) and (2) were used in th’S study. Comparing the results for these two 
viscosities showed that as the rain rate and the turbulent mixing increased, an increase 
in the attenuation of the wave also occurred. The time interval required to dampen a 3 
cm amplitude wave with 10 cm mi.xing depth was 380 wave periods for light rain and 80 
wave periods for heasy rain. Doubling this wave amplitude and maintaining a constant 
10 cm mixing depth, is found to increase the wave periods required to dampen the wave. 
During a light rain event the wave is dampened within 655 wave periods and 1 26 wave 
periods during a hea\y rain. 

b. 20 cm Mixing Depth 

As shown in Table 5, doubling the mixing depth did not double the rate of 
energy change as hypothesized, although a significant (50-80%) change did occur. For 
a wave of 2.8 m wavelength and 6 cm amplitude, a light rain event caused the wave to 
be dampened within 420 wave periods and .vithin 80 wave periods for a hea\y rain event. 
This is a reduction of approximately 38®/o of the wave periods required to dampen the 
same wave with a 10 cm mi.xing depth. 



IV. CONCLUSION 



Strong rain-induced mixing within a ver>' thin surface layer will increase the atten- 
uation rate of surface gravity waves by several orders of magnitude over the rate of at- 
tenuation due to molecular viscosity alone. In this case study a single wavelength (2.8m) 
was used together with two mixed layer depths (10 and 20 cm). In the case of strong 
mixing (v = lO'^m^s-') with a 10 cm thick mixed layer, the wave was dampened with 
12-24 wave periods compared to 46,000 wave periods for molecular viscosity alone. 
Using estimates of the eddy mixing coefficients for hea\7 and light rain from Manton 
(v = 1.6 a- 10-^ s"' and v = S.rlO'^m^s"') , the wave is still dampened within 126 and 

655 wave periods, respectively. Increasing the mixed layer depth to 20 cm produced a 
significant increase in attenuation, reducing the nunber of wave periods required to 
dampen the wave to 80 and 420 wave periods, respectively for hea\w and and light rain. 
This influence of rain on wave attenuation has sev eral potential applications. 

Scatterometer and SAR observations of the sea surface depend on the presence of 
short gravity and capillarv- waves. These waves are generally shorter than the case study 
wavelength (2.8 m) and they will be even more strongly attenuated by rain-induad 
mixed layers. Thus, this theorv' can be used to explain the observed SAR imagerx' in 
rainy regions (dark image) and may allow a correction to be applied to scatterometer 
algorithms used to measure wind speed from small scale ocean roughness. This should 
increase the accuracy of weather forecasts for marine shipping. The damping of slightly 
longer gravity waves present in ship wakes is also of interest. Rain will attenuate these 
waves and thus mask the ship's wake, decreasing the chance of detection by other re- 
mote sensing techiques. Rainfall measurements can also be incorporated into wave and 
surf forecasts as wave length dependent attenuation. This should increase the accuracy 
of such models by decreasing the C;.ergy contained i.; the high frequency spectrum. 

This study has developed a boundary layer that resembles the shape of the water 
surface which may assist in simulating mixing problems more realistically. The eflect 
of turbulent kinetic energy could be introduced numerically by using the turbulent 
kinetic energy equation. This equation would provide a "Reynolds' eddy viscosity coef- 
ficient" for each cell which could then be used directly rather than by defining an arbi- 
trary mixed layer. This will allow variable surface gravity wave attenuation as the mixed 
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layer deepens. Through this method the efTect of a finite duration rain squall on mixed 
layer depth and wave attenuation might be estimated. 

This study focused on one wavelength, two amplitudes and two mixing depths. It 
is recommended that further numerical studies be conducted with larger wavelengths, 
higher wave amplitudes and deeper mixing depths. These factors will identify cutoff re- 
gions where rain-induced attenuation no longer has an effect on the wave. Wave tank 
experiments would also benefit numerical modeling to assist in parameter formations for 
the many inputs necessaiy for turbulent kinetic energy models. 
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